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ABSTRACT 

Previously, Semianalytical Satellite Theories based upon the Generalized Method of Averaging have 
been developed for 

- perturbations with no explicit dependence on time, and 

- perturbations with a strong explicit dependence on time 

While the assumption of time independence (TI) is exact only for zonal harmonics and for static 
atmosphere density models, the assumption has also been applied successfully to develop the 
averaged equations of motion for lunar-solar perturbations of satellite orbits with periods up to two 
days (see AIAA preprints 78-1382 and 75-9). However, recent testing of the lunar-solar short 
periodics produced via the TI assumption for the GPS orbital flight regime ( 1 2 hr period) indicates 
that the relative accuracy of these short-periodics is significantly less than the accuracy of the zonal 
short-periodic variations. 

This paper describes the modifications of the Semianalytical Satellite Theory required to include 
these ‘weak’ time - dependent perturbations. The new formulation results in additional terms in 
the short-periodic variations but does not change the averaged equations of motion. Thus the 
m-monthly terms are still included in the averaged equations of motion. This contrasts with the 
usual approach for the strongly time-dependent perturbations in which the m-monthly (or m-daily, 
if tesseral harmonics are being considered) terms would be eliminated from the averaged equations 
of motion and included in the short-periodics computation. 

Numerical test results for the GPS case obtained with a numerical averaging implementation of the 
new theory demonstrate the accuracy improvement. 


*Section Leader, Space Systems Analysis, Air Force Programs Department 
t Staff Engineer 

**CPT, ARMOR, U.S. Army 
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A SEMI ANALYTICAL SATELLITE THEORY 
FOR WEAK TIME-DEPENDENT PERTURBATIONS 

Outline 

• Review of Analytical Results Tor Time-Independent (TI ) Case 

• Numerical Results for Low Altitude Case w/TI Theory 

t Numerical Results for High Altitude (GPS) Case w/TI Theory 
t Analytical Development of Weak Time-Dependent (WTD) Theory 

• Numerical Results for High Altitude Case w/WTD Theory 
(Zonals, Lunar-Solar, and Solar Pressure) 

• Numerical Results for High Altitude Case w/WTD Theory 
(Zonals, Lunar-Solar, Solar Pressure, and 2x2 Tesserals) 
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SHORT PER1QDICS 
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- SHORT PERIODIC COEFFICIENTS ARE FUNCTIONS OF tHE FIVE SLOWLY VARYING MEAN ELEMENTS 
AND THEREFORE SHOULD ALSO BE SLOWLY VARUNG. 

- COUPLING OF THE FAST VARIABLE SHORT PERIODIC VARIATION WITH THE SEMIMAUOR AXIS 
SHORT PERIODIC VARIATION. 

- FOR CONSERVATIVE FORCES, ANALYTICAL EXPRESSIONS ARE POSSIBLE FOP. eC, AND eD. . 
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LOW ALTITUDE TEST CASE 


• EPOCH CONDITIONS: 1974, Oct. 21,. 

OSCULATING ELEMENTS 

a - 6644.586 
e = .01 

i = 67.98538419° 
ft = 91.99738418° 
w = 200.6741688° 

M = 164.3173126° 

t S/C 
C 0 = 2.0 

? 

Area = 1 .86m 
Mass = 677 , kg 


10 hrs, 24 min. 

MEAN ELEMENTS (PCE) 

a - 6636.3797 
e - .0106045 
T = 67.97090021° 
a = 91.9949106° 
u = 200.21097331° 

M = 164.77124281° 

• ATMOSPHERE 

Modified Harris-Priester 
w/Fio y = 150 


I FORCE MODELS 

COWELL (30 second step) 

< ^ 2 > ****^6 
and drag 


SEMIANALYTICAL (1 day step) 

First Order: an( * Dr *5 

2 

Second Order: + O^-Drag Coupling 

in the AOG 

(IZSAK + Analytical Drag - J ? ) 
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HIGH ALTITUDE TEST CASE 


1. ASSUME A SET OF EPOCH MEAN ELEMENTS; THESE ARE 'CONSTANTS' FOR THE 
SEMI-ANALYTICAL THEORY 

2. AT EPOCH, USE THE SHORT-PERIODIC GENERATOR TO PRODUCE OSCULATING 
ELEMENTS 

3. CONVERT THE OSCULATING ELEMENTS TO POSITION AND VELOCITY; THESE 
ARE THE CONSTANTS FOR THE COWELL THEORY 

4. PROPAGATE THE ORBIT USING .BOTH THE SEMI-ANALYTICAL THEORY AND COWELL 
AND COMPARE THE RESULTING POSITION 'AND VELOCITY HISTORIES 


TEST CASE 9 2 FORCE MODELS 


COWELL 
J 2 J 6 

LUNAR-SOLAR 

SOLAR RADIATION PRESSURE 


SEMI-ANALYTICAL 
J 2 PLUS J 2 2 

LUNAR-SOLAR (TI) 

SOLAR RADIATION PRESSURE (TI) 
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SEMI-ANALYTICAL THEORY 
FOR WEAKLY TIME- DEPENDENT 
PERTURBATIONS 


• OSCULATING EQUATIONS 

da. . 

sr = e F i (a ’ Xtt) 

$ ■ n + cF 6 (J,X,t) 

• ASSUMED FORMS 

da. ± * _ 

~ = e Aj(a,t) a i = a 1 + e n 1 {a , A , t) 

= n^) + e Ag(iT,t) X = X + e T)g(a,\,t) 


• MATCHING EXPRESSIONS FOR da^dt AND dX/dt GIVES 
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THE M-MONTHLIES OUT OF. THE SHORT PERIODICS 

THEN 

i r 2n 




A 1 

H I 

F^(a,X,t) dx, 1 = 

1 


72 



WTD SHORT-PERIODICS 


• DEFINE 


r, S (a,A,t) 


• ASSUME 
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• SUBSTITUTING INTO THE MATCHING EXPRESSIONS GIVES PDE's 
1 3N Io 
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• ASSUME SOLUTION TO PDE 
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• FIRST ORDER RESULT 
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TEST CASE n 


• FORCE MODEL 
COWELL 



LUNAR- SOLAR 

SOLAR RADIATION PRESSURE 


SEMI-ANALYTICAL 

PLUS 0 2 2 

LUNAR-SOLAR (WTD) 

SOLAR RADIATION PRESSURE (WTD) 


• MEAN ELEMENTS 


a = 

26559.5 km 

ft 

= 

o 

O 

o 

e " 

o 

o 

O) 


0.0° 

i = 

O 

w 

o 

0 

M 

* 

o 

o 

o 

• OSCULATING 

ELEMENTS 


• 


a = 

26561.56567 km 

ft 

* 

359.9999657° 

e = 

.00104842 

(0 

= 

359.8560915° 

l = 

63.001124° 

M 

- 

.1436848842° 




TEST CASE U RESULTS 



TIME (DAYS) 

Ax(m) 

Ay(m) 

Az(m) 

RSS(m) 

2 

-.01 

.8 

1.317 

1.54 

4 

-.22 

1.574 

2.704 

3.14 

6 

-.63 

2.278 

4.395 

4.99 

8 

-1.31 

2.866 

6.274 

7.02 

10 

-2.14 

3.601 

7.801 

8.85 

12 

-3.70 

4.702 

9.342 

11.09 

14 

-5.58 

5.322 

10.510 

13.04 
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TEST CASE #3 


» FORCE MODEL 

COWELL 

J 2 ,e<, ”°6 

LUNAR-SOLAR 

SOLAR RADIATION PRESSURE 

(0*5)2 J + ( C » S } g 2 

t MEAN ELEMENTS 

a = 26559.5 km 
e = .001 

i * 63.0° 

• OSCULATING ELEMENTS 

a = 26561.54781 km 
e = .00104802 

1 = 63.001118° 


SEMI -ANALYTICAL 

PLUS J 2 2 

LUNAR-SOLAR (WTD) 

SOLAR RADIATION PRESSURE (WTD) 

(c,s) 2fl + (c,s ) 2>2 


n = o,o° 

w = 0.0° 

M =• 0.0° 


ft = 359.9999706° 
w * 359.8538535° 
M = .1459308175° 
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Figure 1. Radial Difference after 23 hours from Eooch/Semi analytic al 
minus Ccwell for the Lev/ Altitude Circular Test Case 
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Figure 2. Cross Track Difference after 23 hours from Eooch /Semiana Ivti ca 1 
minus Cowell for the Low Altitude Circular Test Case 
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Same initial conditions and perturbations as in Figure i 
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Figure 3. Along Track Difference after 23 hours fron Fpcch/Semina l.yti ca 1 
fnirius Cowell for the Low Altitude Circular Test Case 
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Same perturbations and initial conditions as in Figure 1, 
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Figure 5 Time History of the Drag Sernimojor Axis Short Periodic Coefficients* for tho Low Altitude 
Circular Test Case 

* The 1A coeff icients, Cj j end 


Same initial condition# and perturbation# 
as in drag option 6 of AAS-79-133 
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Figure 6, Radial Difference (n Theory) 
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Figure 7. Cross Track Difference (TI Theory) 


c 

R 

0 

5 

s 

r 

R 

A 

c 

K 

D 

I 

F 

F 

E 

R 

6 
H 
C 
E 
S 

I 

M 

H 

e 

T 

E 

R 

S 


***** POSITION DIFFERENCE 

I. — ■ 


SATELLITE GPS 
LUNAR ONLY (TI) 


11I11U 
.i 


i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i * 
i 

i * 

i* 

i 

i 

i 

i 

i 

i 

i 

> . 
i 
i 
i 
i 

s . 

I 

I 

I 

I 

3 . 

1 

X 

I 

I 


******** 

*** ««* 


*** 

*** 


*** ** 
***** 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

I 

X 

I 

I 

I 

I 

I 

I 

I 

I 

I 


— .1 
15.00 


I.— 

6.0 


1 ' 50 3.00 A. 50 5.06 7.50 9.00 10.50 12.00 13.50 

TIME FROM YYMMDO HHMMSS IN HRS 
790761 0 


0RB1 FILE ON UNIT 2A» DATA RECORDS START AT 790701 0 

QRB1 FILE ON UNIT 51 » DATA RECORDS START AT 790701 0 


USER'S NOTES. 


82 



Figure 8. 


Along Track Difference (TI Theory) 
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Figure 9. Along Track Difference (WTD Theory) 
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Figure 10 o Along Track Difference (WTD Theory) 
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Figure 11* Radial Difference/Semiana lytica 1 minus Cowell for the GPS Test Case 
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Figure 12. Cross Track Difference/Semi'analytical minus Ccwell for the 
GPS Test Case 
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Figure 13. A1 ono Track Difference /Semir.alyti cal minus Cowell for the PPS 
Test Case 
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Same initial conditions and perturbations as in Figure 11. 
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Figure 14. Radial Difference/Semiana lytical -minus Cowell for 
GPS (Test Case #3) 
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Figure 15. Cross Tnack Difference/Semiancllytical minus Cowell for 
GPS (Test Case #3) 
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Figure 16„ Along Track Difference /Semi analytical Sminus Co/ell for 
GPS {Test Case #3) 
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